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Abstract 

A new method for computing all elements of the lattice quark propagator is pro- 
posed. The method combines the spectral decomposition of the propagator, com- 
puting the lowest eigenmodes exactly, with noisy estimators which are 'diluted', 
i.e. taken to have support only on a subset of time, space, spin or colour. We find 
that the errors are dramatically reduced compared to traditional noisy estimator 
techniques. 
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1 Introduction 

Hadron spectroscopy has traditionally been performed in lattice QCD by com- 
puting quark propagators from one or a few points on the lattice (usually the 
origin) to all other points — so-called point or point-to-all propagators — by 
inverting the fermion matrix with a point source, and combining the result- 
ing propagators with appropriate operators to produce the desired hadronic 
correlators. This method does not require massive computing power, but re- 
stricts the accessible physics to primarily the flavour non-singlet spectrum. 
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Flavour singlet mesons, as well as condensates and other quantities contain- 
ing quark loops, require propagators with sources everywhere in space, and 
other methods, notably noisy sources, have been used for this purpose. 

Furthermore, point propagators restrict the interpolating operator basis used, 
for example, to produce early plateaux in effective masses, since a new inver- 
sion must be performed for every operator that is not restricted to a single 
lattice point. 

Point propagators also throw away a large portion of the information contained 
in the gauge configurations. It is likely that a limited number of expensive con- 
figurations with light dynamical fermions will be available in the near future, 
and it will be highly desirable to extract as much information as possible from 
these lattices. 

All-to-all propagators [1,2,3,4,5,6,7,8,9,10,11,12,13] provide a solution to these 
problems, but are usually too expensive to compute exactly as this requires an 
unrealistic number of quark inversions. Stochastic estimates tend to be very 
noisy and variance reduction techniques are crucial in order to separate the 
signal from the noise. In this paper we propose an exact algorithm to compute 
the all-to-all propagator utilising the idea of low-mode dominance corrected 
by a stochastic estimator which yields the exact all-to-all propagator in a 
finite number of quark inversions. Some preliminary results were presented in 



The structure of the paper is as follows: In Sec. 2 we describe the principles 
of the method, and the notation we will use. Sec. 3 describes how QCD cor- 
relators, and meson two-point functions in particular, are computed in this 
framework. In Sec. 4 we present results for the meson spectrum, including 
P-waves and static-light mesons, using this method, and compare these with 
results using traditional methods. Finally, in Sec. 5 we present our conclusions. 



2 Theory and notation 

In the following, latin indices . . . denote colour; greek indices a, /3, . . . de- 
note spin. Indices in brackets denote eigenvectors and dilution indices, while 
indices in square brackets denote independent noise vectors; all of these will 
be introduced in the following. We use (u, to denote the inner product (dot 
product) of two fermion vectors u and v, 



Ref. [14]. 
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If the t-argument is absent from either of the vectors, the product is taken to 
be a global dot product, i.e. summed over all t. 

2. 1 Noisy Estimators and Dilution 

The standard method of estimating the all-to-all quark propagator is by sam- 
pling the vector space stochastically. One generates an ensemble of random, 
independent noise vectors, {?7[i], • • • ,V[N r ]}, with the property 

({V[r](x)<B>V[r](yV))=S x , y , (2) 

where ((• • •)) is the expectation value over the distribution of noise vectors. 
Each component of the noise vectors has modulus 1, 

ri ia (x)*ri ia (x) = 1 (no summation) . (3) 

The solution vectors ip[ r ] are obtained in the usual way, 

i/> [r] (x) =M"V](2/)- (4) 
The quark propagator from any point x to any other point y is given by 

1 N r 

m-Hv,x)% = u [r] ml r] ))% = ^ w E^(yH^y ■ (5) 

This method is noisy because it relies on delicate cancellations in the 0(1) 
noise over many samples to find the signal, which falls off exponentially with 
the separation. We propose to remove the Oil) random noise by "diluting" 
the noise vector in some set of variables (j) such that rj = J2jV^\ resulting 
in an substantial reduction in the variance. A particularly important example 
of dilution for measuring temporal correlations in hadronic quantities is "time 
dilution" where the noise vector is broken up into pieces which only have 
support on a single timeslice each, 

Nt—l 

V(x,t)= ]T Tp\x,t), (6) 

j=0 

where r]^\x, t) = unless t — j. 

Each diluted source is inverted, yielding pairs of vectors, {ip^\r)^}, which 
then gives an unbiased estimator of M' 1 with a single noise source, 

N d -1 

£ ^(x,t)®r]V(x ,t y. (7) 

i=0 
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Fig. 1. The pseudoscalar propagator computed with and without time dilution. 

We show the effect of time dilution on a pseudoscalar propagator on a 12 3 x 24 
lattice in Fig. 1. The circles are the average of 24 noise sources without any 
dilution and the diamonds are from a single noise source which has been time- 
diluted. This particular scenario is analogous to the "wall source on every 
timeslice" method used by the authors of Ref. [15] to estimate the disconnected 
diagrams appearing in hadronic scattering length calculations. Our method 
is, however, more general and can be extended to the spin, colour and space 
components of the source vector. 1 The "homeopathic" limit of the dilution 
procedure, where we have one noise vector for each time, space, colour and spin 
component, results in the exact all-to-all propagator in a finite number of steps, 
because of the property of Eq. (3), see Fig. 2. This limit cannot be reached 
in practice on realistic lattices, but the path of dilution may be optimised 
so that the noise from the gauge fields dominate the errors in the hadronic 
quantities of interest with only a small, manageable number of fermion matrix 
inversions. 



2.2 Spectral Decomposition 



In the following we will make use of the hermitian Dirac operator Q = 75M, 
where M is the usual Dirac operator. The quark propagator is M _1 = Q _1 75, 



A dilution scheme that includes spin and colour, but not time dilution, was pre- 
viously used in Ref. [8]. 
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Fig. 2. A cartoon of possible deviations of the stochastic estimates of the exact 
solution (at Ndu = -/V max ) for different dilution paths. Simply adding noise vectors 
will give a l/y/N behaviour. We have found that simple dilutions typically follow 
the behaviour exhibited by the bottom curve. 

where Q~ l is given by a sum over all eigenmodes, 

N 1 

Q-\ y , x ) = Y,-v (l \y)®v { H^- (8) 

i A i 

Here, Qv® = XiV^' and N is the rank of the matrix. Computing all these 
eigenvectors is clearly not feasible; however, theoretical arguments, supported 
by numerical evidence [16,9,10], suggest that much of the important infrared 
physics in hadronic interactions is encoded in the low-lying eigenmodes. Trun- 
cating the sum in Eq. (8) at some finite N — > N cv should thus yield a good 
estimate of the quark propagator for practical purposes. Such a truncation 
does, however, violate reflection positivity, making it mandatory to correct it. 



2.3 Hybrid Method 



Given the preceding discussion, a natural suggestion would be to try to calcu- 
late as many as possible of the low modes exactly and correct the truncation 
with the noisy method. This gives rise to two concerns: firstly, the correction 
should leave the exactly solved low- lying modes intact; and secondly, it should 
not introduce large uncertanties in the process. We propose that the stochastic 
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method with noise dilution is a natural way to accommodate both of those 
concerns. 



First, we note that the exact N ev low eigenmodes obtained separately naturally 
divide Q into two subspaces, Q — Qo + Qi, defined by 

Qo = E A ^ (i) ®^ )f > Qi= E A^'W j)t . (9) 

i=l j=Nev+l 

Similarly, the quark propagator is broken up into two pieces, Q^ 1 — Q + Q ± , 
where Q Q is the truncated version of Eq. (8) and Q 1 = Q _1 Vi, where V\ is 
the projection operator 

JVcv 

r 1 = i-r = i-j2 v U) <g> v {j)] . (io) 

3=1 

We correct the truncation and estimate Q 1 using the stochastic method, Q 1 = 
((0[r]®?7[ r ])) with N r noise vectors, {f][i], • ■ • ,V[N r ]}- The solutions are given by 

VV] = QlV[r] = Q' 1 (^1%]) • (ii) 

We now apply the idea of dilution to the stochastic estimation of Q l . Each 
random noise vector, r/^, that is generated will be diluted and orthogonalised 
(with respect to v^) so that it can be used to obtain ip[ r y In other words, we 
now have the following set of noise vectors: 

Un {1) ••• n (1) ) ■■■ (n (Nd) ••• n (Nd) )\ 

where the upper indices denote the dilution and the lower indices label the 
different noise samples. Note that the noise vectors are mutually orthogonal 
due to the dilution before an average over different random vectors are taken, 
i.e., 

V[r] ® V$ W, = for alii ^ j. (12) 
This results in smaller variance than the standard method which mixes noise, 
as Eq. (2) shows. 

There is a natural way to combine the two methods to estimate the all-to- 
all quark propagator. The similarity in the structure of Eq. (8) and Eq. (5) 
suggests that one construct the following "hybrid list" for the source and 
solution vectors: 

f,,(l) ,,0Nev) ^1 

w® = \^,---,?T—,VivV,---,ViV™\ (13) 
I A l A NCV J 

u^ = L^\--- ,v^\^\--- ,^ N A (14) 
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Fig. 3. The pion effective mass from 50, 100 eigenvectors and from the hybrid method 
with 100 and a time-diluted noise vector. 

where the indices run over Nhl — N ev + elements. The unbiased, variance 
reduced estimate of the all-to-all quark propagator (for a single random noise 
vector) is then given by 



Nhl 



u®(x,x 4 ) 



(15) 



Note that defining without the projector V\ on the 77-vectors in Eq. (13) 
also gives an unbiased estimator, which may however have a different variance. 
Using the pion as an example, we can demonstrate how positivity is recovered 
from the truncated propagator. In Fig. 3, we show the effective mass from 
the truncated propagator and from the hybrid method with a time, spin, 
colour and space (even-odd) diluted noise vector. The truncated propagator, 
yielding an effective mass that approaches the asymptotic value from below, 
is corrected by the addition of the diluted noisy propagator. 



3 Implementation for QCD correlators 



Besides expanding the range of applications accessible in lattice QCD, all- 
to-all propagators make the construction of hadronic interpolating operators 
considerably simpler. With the (local) meson creation operator 



O r {x,t) = i)(x,t)r<4>(x,t) 



(16) 
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a meson propagator is constructed using traditional, point-to-all propagators 
as follows, 



C r (f,*;0,0) = TrM- 1 (f,t;0,0)r 75 M (x, t; 0, 0) t 7 5 r . (17) 

This construction, where the local operators are traded for non-local objects 
constructed from the quark propagators, may work well as long as we restrict 
ourselves to operators localised at a single lattice site. However, an extended 
source operator involves connecting quark propagators from sites in the vicin- 
ity of the source site. This requires at least an additional inversion and entan- 
gles the matrix inversion and operator construction stages: to a certain extent, 
the matrix inversion requires knowledge of the operators to be used. 

All-to-all propagators eliminate this complication as both source and sink 
operators are constructed purely from local vectors, 

0^ j \x,t)=w[p,t)^ 5 ru$(x,t). (18) 

The extra factor of 75 comes from the use of the hermitian Dirac matrix, 
75M. Two pairs of vectors, (iU[i], W[ij) and (w[ 2 ], U[2]), are needed for the two 
independent noisy estimators of the two propagators. T can now denote an 
extended operator acting on the vector u. For complicated operators such 
as those used to project out hybrid or exotic states, this is a much needed 
simplification. For example, an interpolating operator for the exotic hybrid 
1~ + can be constructed from combinations of gluonic paths projecting out the 
relevant quantum numbers [17]. One term in such a sum may be 

w [i] (xyU z (x)U y (x + e z )U\{x + e y )u [j \x + e y ) , (19) 

where the common t-index is suppressed. A standard P-wave state may for 
example be constructed using 

0^ p (x,t) = w^(x,t)^ 5 (D k u^)(x,t) , (20) 

where 

D kU ( j \x,t) = U k (x,t)u^(x + e k ,t)-u^(x,t) . (21) 

Using all-to-all propagators, correlation functions are also constructed in an 
intuitive manner. Hadronic correlation functions are obtained from interpo- 
lating operators sitting at different time slices, e.g. 

C AB (5t) =EE 0™ A (t)0™ B (t + St) , (22) 



where 



0^ A (t) =^^(x,ty l5 T A u^(x,t) = (w{p), l5 r A u{p)) (23) 
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and 



O&f (t + St) = (w$(t + 5t), l5 T B u^(t + 5t)), (24) 

for isovector two-point correlators. The correlator is constructed as a product 
of these meson operators on different time slices. The operators Y A,B may 
now also include phases for momentum projections. As a consequence, once 
the hybrid list vectors have been constructed, correlation functions for any 
operator can be computed without any additional inversions. 

3.1 Noise recycling 

We can take advantage of having saved different random source/solution sam- 
ples by reusing them in the contraction. In other words, one can generate Nr 
samples of noise vectors, r)[ r ] , save the corresponding solutions, ip[ r ] to disk and 
perform the contraction, 

c AB (t,t ) = ^(«; H (t),r A M[s] (t))( W[s] (to),r B M[r] (t )), (25) 

r<s 

yielding ~ N R samples of the correlation function. The errors correspondingly 
decrease faster than the naive 1/y/Nn, although the measurements are some- 
what correlated. We have seen in our preliminary tests that the error reduction 
is comparable to some dilution choices. It is clear that if one can afford to save 
the noise and solution vectors onto disk, then this is a straightforward method 
of variance reduction for mass-degenerate mesons. 

3.2 Improving Performance 

When constructing one of these hadron operators, such as Eq. (23), on a 
particular time slice t and for a particular operator A, the hybrid technique 
can be manipulated for computational efficency. 

Firstly, compute (and store) 

z$(x,t)= l5 T A u^(x,t),\/j. (26) 

What now remains for the calculation of the time-slice elements of (Eq. 23) is 
N HL x N HL dot-products of spinors on a time slice, since 

0^ A (t) = J2w§(x,ty z $(x,t) ee (w\p),z$(t)) . (27) 

X 
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We note that time dilution is taken as a minimum requirement in the case of 
meson two-point functions and that the hybrid list is naturally divided 
between our two vector spaces, since 



If we evaluate (and store) 




(28) 



(w\p),z$(t)) = ±(v§(t),z$(t)) 



, Vi < iV cv , Vj , 



(29) 



the elements that remain to be computed are 



(^j(t), (1 - P )$(t)) = (v$(t),z$(t)) - £ <*${v$, $(t)) , (30) 



time slice. The second term is simply a weighted sum over terms from the 
eigenvector space of that have already been calculated. The weights are 
given by global dot products of quark vectors which can be calculated exter- 
nally to any particular time slice or operator: Since each 77 has support on 
only a single time slice, the global dot product is equal to the dot product of 
v with the time-slice restricted rj(t). 

This method also reduces the amount of storage required for the noise quark- 
fields by a factor of N t since only the supported time slice needs to be stored 
on disk, all others being identically zero. This argument is easily extended to 
include the conjugate operator and all required momenta. 

Once the general routine to construct a operator field 0^ A (t), as in Eq. (23), 
and the summation over the hybrid list indices required in Eq. (22) are in 
place, the method becomes a black box to the end user. From Eq. (26), we see 
that this user need only supply the subroutines to create the required hadron 
operators from the quark, antiquark and gluon fields on a time slice for each of 
the required operator fields. These routines would only contain the operations 
to be performed on the quark field, such as Djtp, and have no reference to the 
hybrid lists. 

For isoscalar mesons, the disconnected part of the propagator must be in- 
cluded, yielding the following contraction, 



cf- 0) (Mo) = ( w g(t), 7 5r M g(t))( w g(to),75r t M g(to)) 

- (^ff/(t), 75 r M g(t))( w g(t ), 75 rt M g(t )). (32) 



where 



The first term o: 
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3.3 Recipe to compute two-point functions 

1. Calculate N ev eigenvectors, Vj, of 75M with eigenvalues Xj. 

2. Calculate solution vectors, of 

75 M^ W = Vi (33) 

where r\ W are diluted noise vectors, satisfying Eqs (2), (3) and (6) and 

■p 1 = l--p = l-J2 v (j) <g> v (m . (34) 
j'=i 

At least two independent noise vectors need to be generated for the two 
quark propagators required in two-point functions. 

3. With the hybrid lists defined as in Eqs (13) and (14), construct the (time- 
slice) vectors z^(t), 

z^(t)=l^ A u^(t),\/j. (35) 

where the operation r A u9)(t) is defined by the end user, and can also 
represent a conjugate operation. 

4. Construct the operator field C^j A (t) 

0^ A (t) = {wlp)^](t)) 1 (36) 

on a time slice t and store it to disk. Repeat for all time slices and for 
both the source and sink operators. 

5. Contract these operator fields to obtain the two-point correlator C AB , 

Nhl 

C AB (St) =EE 0^ A (t)0^ B (t + St) . (37) 

t i,j 

6. If increased accuracy is required, calculate N cv ' — N cv additional eigenvec- 
tors of 75M. The solution vectors ip® can be projected into the reduced 
orthogonal eigenvector space in 75M by the addition of a correction term, 

j=N cv A 3 

One can also increase the dilution level, from to N' d diluted noise 
vectors, with only N' d —Nd extra quark inversions. For example, if 'j^Mifj = 
V\r\ with r] = 77W + rf 2 \ we need only calculate the solution to 

75 M^ (1) = V X 4 X) (39) 

since — r\ — rf 1 ^ and ip^ — 4> ~ • 
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Fig. 4. Effective masses for isovector mesons with 100 eigenvectors, time, colour, 
spin and space-even-odd dilution and 75 configurations. 

4 Results 



In this exploratory study, we use a data set consisting of 75 quenched config- 
urations at j3 = 5.7 on a 12 3 x 24 lattice. We have used Wilson fermions with 
hopping parameter k = 0.1675, corresponding to m v /m p = 0.50 [18]. The 
100 lowest eigenvectors have been computed for each configuration, and noise 
vectors have been diluted in time, spin, colour and in space on an even/odd 
basis. 



4-1 Isovector mesons 



We have computed correlators for the pion, rho, the l + ~,0 ++ , 1 ++ and 2 ++ 
P-wave states, using spatially extended operators for the P-waves [17]. The 
quark fields were smeared with 10 levels of Jacobi smearing iterations. In 
cases where we use multiple operators for a single state, we have found some 
differences in the overlap depending on the spatial size of the operators where 
typically larger operators were favoured. 



Figure 4 shows effective masses for the pion, rho and P-wave states, using 100 
eigenvectors, time, colour, spin and space-even-odd dilution and 75 configura- 
tions. A lower statistics result for the hybrid 1~ + was reported in Ref. [14]. 
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Fig. 5. Comparison of spectrum from all-to-all and point-to-all propagators. The 
open symbols are point-to-all, the closed symbols all-to-all with 100 eigenvectors 
and time, colour, spin and space-even-odd dilution. 

We can also compare the results using our method with results using tradi- 
tional point-to-all propagators on the same number of configurations. These 
are shown in Fig. 5 for the pion, rho and l +_ P-wave. The P-wave was com- 
puted using point sources at 3 different spatial points in order to construct 
the derivative operator. While the improved signal for the pion may not jus- 
tify the additional computational expense, it is clear that for the rho meson 
and in particular the P-wave the reduction in noise is very significant with 
all-to-all propagators. This is particularly important if only a limited number 
of configurations is available, which will be the case for light, dynamical full 
QCD configurations. In the case of P-waves (as well as hybrids, and of course 
isoscalar mesons) all-to-all propagators may be the only way of obtaining any 
acceptable signal at all. 



4-1-1 Effect of dilution 

First we study the effect of dilution using only noise vectors and no eigenvec- 
tors. Results for the pion, rho, 1 ++ and 2 ++ are shown in Figures 6-10 
with square symbols. The error bars on the errors were obtained by bootstrap 
resampling. Time-dilution is sufficient to reach the level where the error on 
the correlator is saturated by the gauge fluctuations for the pion. This is not 
the case for the vector as the error decreases with increasing number of dilu- 
tions. We also show the curve which would be expected to result from using 
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Fig. 6. Demonstration of the effect of dilution on the statistical errors as a function 
of dilution, for the pion. The squares denote the errors obtained using only noise 
vectors, while the diamonds are errors obtained using the hybrid method. 

N independent noise samples per time slice. This curve is given by 

9 0~\ — CT? 9 , . 

<>% = - ± J^ JL + ° 9 , (40) 

where a g is the gauge noise, i.e. the statistical uncertainty resulting from 
having a finite number of gauge configurations. We have assumed that the 
gauge noise is given by the errors from the hybrid method with our highest 
dilution level, a g = a^f- We see in all cases studied here that dilution always 
works better than accumulating an equivalent number of independent noise 
samples. A closer inspection reveals that the errors fall off like 1/N rather 
than the naive 1 / This is in addition to the exponential gain due to time 
dilution. 



4-1-2 Effect of eigenvectors 

Figures 3 and 11 shows how the low- lying eigenvectors saturate the signal for 
the pion and to a lesser extent the rho meson. For early times the effective 
mass from the eigenvectors alone clearly exhibits violation of positivity, but 
from t = 7 onwards for the pion and t = 9 for the rho, the full signal for the 
effective mass is captured by the eigenvectors. We also note the clear difference 
between the data from 50 and 100 eigenvectors. We have not performed any 
systematic study of the effect of varying the number of eigenvectors in the 
hybrid method, since this is expected to be highly dependent on the action 
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Fig. 7. Demonstration of the effect of dilution on the statistical errors as a function 
of dilution, for the p. The squares denote the relative errors on the correlator on time 
slice 3 obtained using only noise vectors, while the diamonds are errors obtained 
using the hybrid method. The solid line is the expected behaviour of the errors from 
accumulating an equivalent number of independent (time-diluted) noise samples. 
The inset shows the relative errors from the hybrid method in more detail. 
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Fig. 8. As Fig. 7, for the 1 H P-wave. 
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Fig. 9. As Fig. 7, for the 1 ++ P-wave. 
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Fig. 10. As Fig. 7, for the 2++ P-wave. 

and lattice spacing used, and the non-chiral action and coarse lattices used 
in this exploratory study would be unlikely to yield valuable information for 
more realistic simulations. 



The diamonds (lower points) in Figs 6-10 show the behaviour of the errors 
for the pion, rho, l + ~, 1 ++ and 2 ++ as a function of the dilution level for the 
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Fig. 11. The p effective masses from the 50 and 100 lowest eigenvectors alone, and 
from the hybrid method with 100 eigenvectors. The corresponding figure for the 
pion is shown in Fig. 3. 

hybrid method with 100 eigenvectors. We observe a dramatic decrease in the 
relative error of the estimated correlator when the lowest 100 eigenvectors are 
used. The exception is the pion, where the errors remain constant at about 
5% (from our 75 configurations) regardless of dilution and eigenvectors. We 
interpret this as an effect of the pathologies of Wilson fermions on coarse 
lattices with light quarks, where the large fluctuations of the lowest eigenmode 
affect the pion in particular. 



The insets of Figs 7 and 8 show the errors for the hybrid method in more 
detail as a function of dilution, for the p meson and 1 H P-wave. The gauge 
noise level appears to be reached somewhere between dilution levels 6 (colour 
and space even/odd) and 12 (colour and spin), although a slight downward 
trend may still be observed at A^u = 24 for the P-waves. It is interesting to 
note that fractional error reached for the vector meson is much smaller than 
the pseudoscalar (for this rough action). 



4-2 Isoscalar mesons 



We have computed the disconnected part of the pseudoscalar isoscalar (rf) 
meson. Figure 12 shows the errors as a function of dilution level for the hybrid 
method (N ev = 100) and using only a noisy estimator. We see that the errors 
from the gauge noise are large, of the order of 20%, and consequently the 
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Fig. 12. As fig. 7, for the disconnected piece of the isoscalar correlation function. 

total uncertainty is saturated by the gauge noise almost immediately. This 
is not surprising on this coarse lattice with a non-chiral action, but it gives 
confidence that a good signal can be obtained from a proper action. 



4-3 Static-light mesons 



The benefit of using all-to-all propagators is particularly evident in simulations 
of heavy-light mesons in the static limit. Using a single point propagator for 
the light quark means that the source and sink of the static-light correlator 
are restricted to a single spatial lattice site. All-to-all propagators allow us to 
place source and sink operators at each spatial site on the lattice, yielding a 
dramatic increase in statistics. The application of all-to-all propagators to the 
static-light spectrum has previously been examined in Refs. [6,19]. 

The mesons considered in the previous sections are classified according to their 
J PC quantum numbers. Here, however, the mesons contain non-degenerate 
quarks and charge conjugation is no longer a symmetry of the meson. Also, 
static-light mesons which differ only in the spin of the static quark are de- 
generate. Therefore, in the static limit, we find a single S-wave channel and 
two distinct P-wave channels. We label these Jf where Jg is the total angular 
momentum of the light degrees of freedom. The S-wave is then | and the 
P-wave channels are | + and | . 

Figure 13 shows effective masses for the static-light mesons. The run parame- 
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Fig. 13. Effective masses for the static-light S-wave and P-waves, calculated using 
100 eigenvectors and time, colour, spin and space-even-odd dilution. 



ters are the same as for the isovector mesons shown in Fig. 4 and, once again, 
we have used extended operators for the P-waves. We obtain excellent signals 
in each channel and, as Fig. 14 shows, we are able to determine the P-wave 
ground-state energies to within 1% accuracy. Physically, the | + is expected 
to be heavier than the | + [19] but, given the coarseness of the lattice and the 
actions used, we do not expect our results to be physically meaningful. 

The results were obtained using the variational method with five operators 
differing by the Jacobi smearing parameters of the light quark. This is a very 
simple way of extending the basis operators to get a better overlap with the 
ground state and for identifying excited states. As Fig. 15 shows, we find clear 
signals for the first and second excited states of the S-wave. We can also resolve 
signals for the first excited states of both P-waves. 

Figures 16-18 show the relative errors for the static-light meson correlators as 
a function of dilution level. For the | and | + the errors are already saturated 
by the gauge noise with time dilution alone, while for the | + there is still some 
gain to be had from additional dilution. In this case we find that the errors fall 
off approximately as 1/x/iV, indicating little if any sensitivity to the dilution 
path for this variable: after time dilution, the purely statistical gain is the 
dominant factor. 
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Fig. 14. The energy-difference between the static-light P-wave ground states, with 
the line denoting the best fit. 
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Fig. 15. Effective masses for the three lowest-lying states of the static-light S-wave. 
5 Conclusions 



We have presented a new algorithm to estimate the all-to-all propagator. All- 
to-all propagators make it possible to make use of all the available information 
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Fig. 16. Demonstration of the effect of dilution on the error-bars as a function of 
dilution, for the static-light S-wave, with 100 eigenvectors. The solid line is the 
expected behaviour of the errors from accumulating an equivalent number of inde- 
pendent, time-diluted noise samples. 
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Fig. 18. As fig. 6, for the static-light | P-wave. 

in a gauge configuration, which considering the cost to generate full QCD con- 
figurations may be of crucial importance. They are also a necessary ingredient 
in flavour-singlet physics and in computing fermionic thermodynamic quanti- 
ties. 

All-to-all propagators have a further advantage over point propagators in that 
operator construction is considerably simplified: the operators are constructed 
in a natural way from local fields, and extended operators used in variational 
methods may be employed at no additional cost. 

We have presented evidence that diluting stochastic estimators in time, colour, 
spin or other variables gives less variance than traditional noisy estimators. 
This is not unexpected, since dilution will yield the exact all-to-all propagator 
in a finite number of iterations. More work is needed to determine the optimal 
dilution path for different observables, although the simplest choices (colour, 
space even/odd, and to some extent spin dilution) are often sufficient to reach 
the gauge noise level for the ensemble sizes considered here. 

The hybrid method allows one to extract the important physics from low-lying 
eigenmodes and combine this with a noisy correction in a natural way. This 
may be implemented so that the user can be blind to the details of the dilution 
and eigenvector list. 

We have found that this method makes it possible to determine traditionally 
"noisy" quantities such as masses of P-wave states, hybrids and static-light 
mesons (including excited states) to a high level of precision with limited 
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statistics. 



Although we have focused on (mainly isovector) mesons in this study, the 
method is also straightforwardly applicable to baryons and to thermodynamic 
quantities (condensates and susceptibilities). Work on flavour singlets and 
thermodynamic observables is currently in progress. 
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